Loschmidt cooling by time reversal of atomic matter waves 
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We propose an experimental scheme which allows to realize approximate time reversal of matter 
waves for ultracold atoms in the regime of quantum chaos. We show that a significant fraction of 
the atoms return back to their original state, being at the same time cooled down by several orders 
of magnitude. We give a theoretical description of this effect supported by extensive numerical 
simulations. The proposed scheme can be implemented with existing experimental setups. 
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The statistical theory of gases developed by Boltzmann 
leads to macroscopic irreversibility and entropy growth 
even if dynamical equations of motion are time reversible. 
This contradiction was pointed out by Loschmidt and is 
now known as the Loschmidt paradox The reply of 
Boltzmann relied on the technical difficulty of velocity 
reversal for material particles [2|: a story tells that he 
simply said "then go and do it" . The modern resolution 
of this famous dispute came with the development of the 
theory of dynamical chaos 0, 0- Indeed, for chaotic 
dynamics small perturbations grow exponentially with 
time, making the motion practically irreversible. This 
explanation is valid for classical dynamics, while the case 
of quantum dynamics requires special consideration. In- 
deed, in the quantum case the exponential growth takes 
place only during the rather short Ehrenfest time [6J] , and 
the quantum evolution remains stable and reversible in 
presence of small perturbations Q. Quantum reversibil- 
ity in presence of various perturbations has been actively 
studied in recent years and is now described through the 
Loschmidt echo (see @ and Refs. therein). This quantity 
measures the effect of perturbations and is characterized 
by the fidelity f(t r ) = \(tp p (2t r )\ip(0))\ 2 , where \tp p ) is 
the time reversed wavefunction in presence of perturba- 
tions, \ip) is the unperturbed one and t r is the moment 
of time reversal. Experimental implementations of time 
reversibility for quantum dynamics or propagating waves 
have been realized with spin systems (spin echo tech- 
nique) 0, acoustic [l(| and electromagnetic 11 1 waves, 
resulting in various technological applications. Surpris- 
ingly enough, the reversibility signal becomes stronger 
and more robust in the case of chaotic ray dynamics [lfj . 
However, despite the significant experimental progress 
made recently in the control of quantum systems, the 
time reversal of matter waves has not been performed so 
far. 

Here we present a concrete experimental proposal of 
an effective time reversal of atomic matter waves. The 
proposal relies on the kicked rotator model, which is a 
cornerstone model of quantum chaos 0, 0, [Hj]. This 
model has been built up exp erimen tally with cold atoms 
in kicked optical lattices [13, [lj, [H, [l6[ • Recent progress 
allowed to implement this model with ultracold atoms 



and Bose-Einstein Condensates (BEC) [1J, [H, Ha l2Q|, 
with high-precision subrecoil definition of the momentum 
of the atoms, allowing for example to observe 17, n| high 
order quantum resonances 2l|. We show that these ex- 
perimental techniques allow to perform time reversal for 
a significant part of the atoms. Surprisingly, this fraction 
of the atoms becomes cooled down by several orders of 
magnitude during the process. We call this new cooling 
mechanism Loschmidt cooling since it is directly related 
to the time reversibility. 

The quantum kicked rotator corresponds to the quan- 
tization of the Chirikov standard map [1, [|| : 



P = P 



x + Tp 



(1) 



where x is the position and p the momentum of an atom, 
and bars denote the variables after one map iteration. 
Here a; is a continuous variable in the interval ] — oo, +oo[. 
The physical process behind corresponds to rapid change 
of momentum created by a kick of optical lattice followed 
by a free propagation of the atoms between periodic kicks 
of period T. The classical dynamics depends only on the 
single parameter K = kT ', and undergoes a transition 
from integrability to chaos when K is increased. Global 
chaos sets in for K > K c = 0.9716.... The dynamics of 
(fT|) is time reversible, e.g. by inverting all velocities at 
the middle of the free motion between two kicks. 

The quantum evolution over one period is described by 
a unitary operator U acting on a wavefunction: 



ip = Uip 



-iTp /2 — ik cos x 



1> 



(2) 



where the momentum p is measured in recoil units, and 
T plays the role of an effective Planck constant. The 
momentum operator p = —id/dx has eigenvalues p = 
n + (3 where n is an integer and f3 is the quasimomentum 
for a wave propagating in the x direction. The particle 
energy is E = E r p 2 /2 where E r is the recoil energy. We 
consider noninteracting atoms. It is convenient to express 
the time t in number of map iterations. In experiments, 
values of time up to t — 150 have been achieved [l6j |. 
Also, a very narrow initial momentum distribution down 
to rms erg w 0.002 can be reached with BEC [l7j |. Values 
as high as k « 4 have been realized experimentally with 
T varying between 1 and 4-7T [l|| . 
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FIG. 1: (Color online) Dependence on time t of the rescaled 
IPR £t/£o (solid curve with blue/black squares) and of 
the rescaled inverse probability at zero momentum Wo/Wt 
(dashed curve with red/gray dots). Initial state is a Gaussian 
wave packet with rms op — 0.002 and kBTo/E r = op/2 = 
2 x 10 -6 (see Fig.|5J. Here k = 4.5, T = 4?r + e for < t < U 
and T' = 4tt - e for t r < t < 2t r with e = 2 and t r = 10. 

To perform the time reversal, we first write T as T = 
4ir + e. After t r iterations of (O, we interchange the order 
of kick and free propagation, change T to T' = 4ir — e 
and k to —k, and let the system evolve during another t r 
iterations. Such a modification of T can be easily realized 
experimentally (see e.g. [l7j| ) . The sign of k can be in- 
verted by changing the sign of the detuning between the 
laser and the atomic transition frequencies or through a 
shift of the optical lattice by half a wavelength. Then 
for t > t r the map © becomes U r = e ik cos * e~ iT ' ^ I ' 2 
where the second operator acts on momentum eigenstates 

\ n + m as e ~iT' (n+Pf/2 _ e tT(n+P) 2 /2 e -8i7r/3(n+/3/2) _ 

Thus U r = C/t e -8"r/3(ft+/3/2) and the components with 
(3 = (integer values of p) are exactly reversed, while 
for other small (3 values the time reversal works only ap- 
proximately. This determines the meaning of approxi- 
mate time reversal. Another way of inverting k is to 
use in the first and last propagation steps for t > t r the 
value T" = 2ir — e instead of T". Indeed for (3 = wc 

. -2 i -.-2 -i 

have e mn e - lKcosx e lm — e"™ 1 , and thus inversion of 
k works again exactly for (3 = and approximately for 
small (3 values. In the following, we use the first method 
of k inversion, but we checked that the second method 
gives similar results. 

To characterize the wave packet dynamics, we use 
the Inverse Participation Ratio (IPR) denned by £t = 
£pWp(t)/E p W?W where W p (t) = \(p\m)\ 2 is the 
probability in momentum space at time t. Another use- 
ful quantity is the probability at zero momentum Wt = 
Wo (t) . For numerical simulations we used up to N = 2 23 
discretized momentum states with Ap = 1/40000. Both 
quantities £ t and W t are rescaled by their values at t = 
that makes them independent of Ap. Their dependence 
on time is shown in Fig. Q] which clearly demonstrates 
time reversal of both quantities. The return for £ t is not 
perfect due to the contribution of nonzero quasimomen- 
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FIG. 2: (Color online) Initial states probability (dashed 
curves) and final return probability (solid curves) shown as 
a function of quasimomentum (3. Here fcsTb/ ' E r = 2 x 10~ 4 
(red/gray curves) and fcsTo/i? r = 2 x 10~ 6 (blue/black 
curves), the latter corresponding to data of Fig. [T] and be- 
ing similar to experimental conditions [l^. For return prob- 
abilities, the central peaks coincide for the two values of To. 
Insets: probability W v at t = t r = 10 (left) and final probabil- 
ity W p at t = tf = 2t r (right) on a larger scale with p = n + (3 
for ksTo/Er = 2 x 10 -6 . All distributions are scaled by the 
value of the initial probability at zero momentum ((3 = 0). 
Parameters are as in Fig. [T] 



turn components, whereas the curve for Wt is perfectly 
symmetric. We note that £t shows a diffusive growth for 
t < 4 followed by a saturation due to dynamical localiza- 
tion @, 0,0, Q for t < t r . 

The momentum probability distributions at initial t = 
and final t = tj = 2t r times are shown in Fig. [5] 
The striking feature is that the final distribution in 
—0.5 < p < 0.5 is much more narrow than the initial one 
and has the same maximal value since the probability at 
(3 = returns exactly. The shape of the reversed peak is 
independent of Tq. This can be interpreted as a cooling of 
the atoms which remained in this momentum range, that 
defines the Loschmidt cooling mechanism. The narrowing 
of the central peak means that in compensation a signif- 
icant fraction of atoms has obtained higher momentum 
values p as is clearly seen in the right inset of Fig. O But 
even if the full distribution is rather broad the reversed 
peak is clearly dominant. On the contrary, the left inset 
showing the distribution at t = t r displays homogeneous 
chaotic distribution of momentum components. Thus it 
is the time reversal which produces the peak at the origin 
and performs effective cooling. It is natural to define the 
size of the return Loschmidt peak by its half width (3l 
with Wp L (2t r ) = Wp=o(2t r )/2. The fraction of returned 

Oft 

atoms is Pp = Y^-2/3 L Wp(2t r ) and their temperature is 
T f = Y?% L (3 2 W0(2t r )/(2P,3). Similarly to the case of 
chaotic acoustic cavities 10 . quantum chaos makes the 
time reversed peak more visible. From an experimental 
viewpoint, the atoms outside the reversed peak can be 
eliminated by an escape process while those inside can 
be kept by switching on a suitable trap potential or at- 
traction between atoms (e.g. Feschbach resonance) . Such 
a procedure is similar to the process of evaporative cool- 
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FIG. 3: (Color online) Density plot of the return probability 
distribution W p (2t r ) as a function of the rescaled atom en- 
ergy E/ksTo and of the kick strength k, where E — E r p 2 /2 
and k B T /E r = 2 x 1CT 6 ; here e = 2 and t r = 10. Colors 
denote density from white (minimal) at the right to red/gray 
(maximal) at the left. Black curve shows the final tempera- 
ture Tf of the return Loschmidt peak (E — > fcgT/ ■ = (E)) as 
a function of k. Logarithm is decimal. 



ing. However in our case the effective evaporation takes 
place very rapidly due to dynamical chaos. 

The variation of the final return probability distribu- 
tion in energy E = E r p 2 jl with k is shown in Fig. [3] 
When k increases the distribution in energy becomes 
more and more narrow, so that Loschmidt cooling be- 
comes more efficient. This is related to the fact that the 
dynamics becomes more chaotic as k increases. The tem- 
perature Tf drops by almost two orders of magnitude, 
showing significant oscillations with k. 

The decrease of Tf/To with k is shown in more de- 
tails in Fig. |U It is related to the increase of the local- 
ization length I of quasienergy eigenstates with k. In- 
deed, it is known that I = D q /2 where D q = k 2 g(K q )/2 
is the quantum diffusion rate, K„ = 2fcsin(T/2) be- 
ing the quantum chaos parameter 0, H3] . The function 
g{K q ) takes into account the effects of quantum corre- 
lations and is given [22| by g{K q ) = for K q < K c , 
g{K q ) « 0A2(K q - K c ) 3 /K 2 for K c < K q < 4.5 and 
g(K q ) = 1 - 2J 2 (K q ) - 2J 2 (K q ) + 2J$(K q ) + 2J|(X g ) 
for K q > 4.5 where J m are Bessel functions. Due to 
this localization there is always a residual probability 
Wp ~ l/lcrp m the interval —0.5 < p < 0.5, even in ab- 
sence of time reversal. However this residual probability 
is much smaller than the maximum of the reversed peak 
Wq ~ I/073. The width of this peak can be estimated 
as follows: (3 > in U r acts as a small perturbation of 
the exactly reversible operator, whose eigenstates have 
M ~ I components. This perturbation gives after time t r 
an accumulated quantum phase shift in the wavefunction 
Acj) = 8n[3nt r «s &Tt(3Mt r . Thus only the atoms with f3 < 
Pl ~ l/(8%Mt r ) return to their initial state, and their 
fraction is Pp ~ Pl/v^- The Loschmidt temperature of 
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FIG. 4: (Color online) Ratio of final to initial temperatures 
(solid curves) as a function of k, from top to bottom for 
k B T /E r = 2 x 10~ 6 (blue/black), k B T /E r = 1.8 x 10~ 5 
(green/light gray), and fcsTo/ ' E r = 2 x 10~ 4 (red/gray), with 



2 and t r 



10. 



Dashed curves show the theory ©. Full 
2 x 10~ 4 , confirming theory 
(|4|). Inset : same ratio as a function of e for k — 4.5; solid 
curve shows numerical data, dashed curve is the theory @. 



circles show Pp for ksTo/E,, 



these atoms is k B T L = S r /3|/2 w E r / '(1287T 5 'CD\t 2 r + 1) 
where C is a numerical constant which according to our 
data is C « 0.4. Thus the ratio T f /T is: 

T f /T = T L /(T + T L ) , k B T L « E r /(500D 2 q t 2 r + 1) (3) 

The formula ([3]) interpolates between the weakly pertur- 
bative regime with I -C 1 and the strong chaos regime 
with I » 1. This theory assumes that k B To -C -Err- 
and that the localization time scale t* « D q is shorter 
than t r which is approximately satisfied for most k val- 
ues in Fig. |4j The theory ([3]) satisfactorily describes 
the global behavior of Tf/To as can be seen in Fig. 2] 
Small scale oscillations should be attributed to meso- 
scopic fluctuations. These fluctuations are stronger when 
e is varied (see inset of Fig. HJ) , that should be attributed 
to high order quantum resonances 2l| at rational val- 
ues of T/47T [23j. We checked that the cooling remains 
robust even in presence of 1% fluctuations of e during 
map iterations. For the case t r <c t* one should use 
M rs yfDqtr instead of M = D q since the diffusion takes 
place during the whole time interval t r . In this case, 
k B T L « E r /{128Tr 2 CD q t 3 r + 1). 

It is important to evaluate the fraction Pp of returning 
atoms. The estimates given above lead to Pp ~ PlI&p oc 
l/(apD q t r ) and thus: 



Tf/To « P, 



0' 



(4) 



that is well confirmed by the data in Fig. 2] The formula 
Q is written for dimension d = 1. For higher dimen- 
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e.g. 0, 0, El, 0, E3, EE EH S3)- The realization of 
our proposal will shed a new light on the long-standing 
Boltzmann-Loschmidt dispute on time reversibility. 

We thank CalMiP for access to their supercomputers 
and the French ANR (project INFOSYSQQ) and the EC 
project EuroSQIP for support. 
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FIG. 5: (Color online) Fidelity / (squares), fraction Pp of 
atoms in the Loschmidt peak (triangles) and fraction P x of 
these atoms in the coordinate space cell (circles) (see text) as 
a function of time reversal moment t r . Straight dashed curve 
shows the fit logPp = —1.13 log t r + 0.43. Parameters are as 
in Fig. [TJ Logarithms are decimal. 



sions it generalizes to Tf/T w P 2 J G 1 . The data displayed 
in Fig. [5] show that Pp drops approximately as l/t r , in 
agreement with the estimate above. For example, the 
cooling by a factor of 100 is reached for 10% of atoms. 
We also checked that a significant fraction of the atoms 
returns not only in momentum space but also in coordi- 
nate space. To this end we define the fraction of atoms 
P x which are in the return Loschmidt peak and in the 
coordinate space interval [— 2/ap, 2/ap]. The data show 
that P x is close to the value of Pp for moderate values of 
t r < 10, so that the Loschmidt cooled atoms remain close 
both in momentum and coordinate space and thus can 
be efficiently captured by a trap potential or a Feschbach 
resonance. We note that at large t r values Pp can become 
smaller than the fidelity / = \(if>(0)\i/j(2t r ))\ 2 , which char- 
acterizes proximity of the whole wavefunctions and not 
only the reversed peak. 

In conclusion, we have presented a concrete experimen- 
tal proposal of time reversal of matter waves of ultracold 
atoms in the regime of quantum chaos. If the atoms 
were classical particles, they would never return back due 
to exponential instability of dynamical chaos. But the 
quantum dynamics is stable and thus a large fraction of 
the atoms returns back even if the time reversal is not 
perfect. This fraction of the atoms exhibits Loschmidt 
cooling which can decrease their temperature by several 
orders of magnitude. The reversed peak is very sensitive 
to variations of (3 and other parameters breaking time 
reversal symmetry, and therefore this setup can be used 
as a sensitive Loschmidt interferometer to explore such a 
symmetry breaking (e.g. a gravitational field component 
along the optical lattice gives a shift in (3 which affects 
the Loschmidt peak). The parameters considered here 
are well accessible to nowadays experimental setups (see 
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